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A central goal in condensed matter and mod- 
ern atomic physics is the exploration of many- 
body quantum phases and the universal charac- 
teristics of quantum phase transitions in so far 
as they differ from those established for thermal 
phase transitions. Compared with condensed- 
matter systems, atomic gases are more precisely 
constructed and also provide the unique opportu- 
nity to explore quantum dynamics far from equi- 
librium. Here we identify a second-order quan- 
tum phase transition in a gaseous spinor Bose- 
Einstein condensate, a quantum fluid in which su- 
perfluidity and magnetism, both associated with 
symmetry breaking, are simultaneously realized. 
8 'Rb spinor condensates were rapidly quenched 
across this transition to a ferromagnetic state and 
probed using in-situ magnetization imaging to 
observe spontaneous symmetry breaking through 
the formation of spin textures, ferromagnetic do- 
mains and domain walls. The observation of topo- 
logical defects produced by this symmetry break- 
ing, identified as polar-core spin- vortices contain- 
ing non-zero spin current but no net mass cur- 
rent, represents the first phase-sensitive in-situ 
detection of vortices in a gaseous superfluid. 

Most ultracold atomic gases consist of atoms with non- 
zero total angular momentum denoted by the quantum 
number F, which is the sum of the total electronic angu- 
lar momentum and nuclear spin. In spinor atomic gases, 
such as F = 1 and F = 2 gases of 23 Na 1, 2] and 87 Rb 
0. [H ■ all magnetic sublevels representing all orientations 
of the atomic spin may be realized @. The phase co- 
herent portion of a Bose-Einscin condensed spinor gas 
is described by a vector order parameter and therefore 
exhibits spontaneous magnetic ordering. Nevertheless, 
considerable freedom remains for the type of ordering 
that can occur. For 87 Rb F = 1 spinor gases, the spin- 
dependent energy per particle in the condensate is the 

sum of two terms, C2n(F} 2 + q(F 2 ), where F denotes the 
dimensionless spin vector operator. The first term de- 
scribes spin-dependent interatomic interactions, with n 
being the number density and C2 = {Aith 2 /3m) (0,2 — do) 
depending on the atomic mass m and the s-wave scat- 
tering lengths a / for collisions between pairs of particles 
with total spin / H, 0- Given C2 < for our system 
0, EL HI 0- t ne interaction term alone favors a ferro- 
magnetic phase with broken rotational symmetry. The 
second term describes a quadratic Zeeman shift in our ex- 
periment, with q = {h x 70Hz/G 2 )£? 2 at a magnetic field 
of magnitude B [l(| ■ This term favors instead a scalar 
phase with no net magnetization, i.e. a condensate in the 



\m z = 0} magnetic sublevel. These phases are divided by 
a second-order quantum phase transition at q — 2\c2\n. 

This article describes our observation of spontaneous 
symmetry breaking in a 87 Rb spinor BEC that is rapidly 
quenched across this quantum phase transition. Nearly- 
pure spinor Bose-Einstein condensates were prepared in 
the scalar \m z = 0) phase at a high quadratic Zeeman 
shift (q ^> 2|c2|n). By rapidly reducing the magnitude 
of the applied magnetic field, we quenched the system 
to conditions in which the ferromagnetic phase is ener- 
getically favored (q -C 2|c2|n). At variable times TLqJj 
after the quench, high-resolution maps of the magnetiza- 
tion vector density were obtained using magnetization- 
sensitive phase contrast imaging Soon after the 
quench, transverse ferromagnetic domains of variable size 
formed spontaneously throughout the condensate, di- 
vided by narrow unmagnetized domain walls separating 
domains of nearly opposite orientation. Concurrent with 
the formation of these domains, we also observed topo- 
logical defects that we characterize as single charged spin 
vortices with non-zero circulation of spin current and an 
unmagnetized filled core. 

Spinor BECs in the \F = l,m 2 = 0) hyperfine state 
were confined in an optical dipole trap characterized 
by oscillation frequencies (u> x , oj y , u> z ) = 27r(56, 350, 4.3) 
s _1 . The condensates, typically containing 2.1(1) x 10 6 
atoms, were formed at a magnetic field of 2 G and were 
characterized by a peak density no — 2.8 x 10 14 cm~ 3 
and Thomas-Fermi radii (r x ,r y ,r z ) = (12.8,2.0,167) 
fj,m (Supplementary Note 1). The anisotropic conden- 
sates were effectively two-dimensional with respect to 
the spatial variations in the internal-state wavefunc- 
ti on due to th e fact that the spin healing length, £ s = 
yj h 2 / 2m\c2\nQ = 2.4 fim, is larger than r y , therfore elim- 
inating spin dynamics in the y direction. Thus, imaging 
the condensate in the x—z plane produces complete maps 
of the magnetization density. The reduced dimensional- 
ity also determines the nature of topological defects that 
may arise in this quantum fluid. 

After the condensate was formed, the magnetic field 
was ramped linearly over 5 ms to a magnitude of 50 mG 
in the i direction and was held at this setting for a vari- 
able time T^qJj prior to the imaging sequence. At this 
field, the quadratic Zeeman energy q — hx 0.2 Hz is negli- 
gible compared to the spin-dependent interaction energy 
of (4/5)2|c2|no = h x 8.2 Hz, which is computed using the 
density averaged in the y direction. 

The condensate magnetization was measured in-situ 
using phase-contrast imaging with pulses of circular po- 
larized light, yielding an optical signal given approxi- 
mately as C (l + l{F y ) + l(F 2 )^j 03. Here, C = ^0-7/8^ 



2 




1R RR 



15B 1RR 1SR 91R 



FIG. 1: Direct imaging of inhomogeneous spontaneous mag- 
netization of a spinor BEC. Transverse imaging sequences 
(first 10 of 24 frames taken) are shown (a) for a single conden- 
sate probed at 2j, ih = 36 ms and (b) for a different conden- 
sate at r n old = ^ ms. Shortly after the quench, the system 
remains in the unmagnetized \m z = 0) state, showing nei- 
ther short-range spatial nor temporal variation (i.e. between 
frames). In contrast, condensates at longer times are spatially 
inhomogeneous and display spontaneous Larmor precession as 
indicated by the cyclical variation of signal strength vs. frame 
number. Orientations of axes and of the magnetic field are 
shown at left. 




where n is the column density of the gas, a = 3A 2 /2tt is 
the resonant cross section, A=795 nm is the wavelength of 
the probe light, 5 is the detuning from resonance, and 7 
is the natural linewidth (Supplementary Note 2). Thus, 
given the column density of the gas and ignoring the 
small (Fy) signal, the phase-contrast signal is an instan- 
taneous measure of one component, F y = (Fy), of the 
(dimensionless) magnetization of the gas. 

We determined all three components of the vector 
magnetization density with repeated images of the same 
atomic sample. First, a sequence of 24 images was taken 
for which the magnetic field remained oriented along the 
z direction (Supplementary Note 2). Larmor precession 
of the transverse magnetization at around 35 kHz is im- 
aged as an aliased low-frequency oscillation (Fig. 1). The 
complex transverse magnetization Ft — F x + iF y can 
then be determined from the relation A(p) exp(z0(p)) = 
«(5/6) Q{p) F T (p) where A(p) and (f>{p) are the amplitude 
and phase of this oscillation at each pixel position p. 
The normalization constant was obtained by aver- 
aging the signal over all image frames and fitting with 
a Thomas-Fermi distribution for the density of the con- 
densate. Following this sequence, two additional image 
frames were obtained within 5 ms in which the magnetic 
field was adiabatically reoriented in the y and —y direc- 
tions. The longitudinal magnetization was determined 
from the difference between these last frames. 

As shown in Fig. 1, at short times after the quench 
(Tjiold < 50 ms), the transverse magnetization images 
show no significant variation across the cloud or between 
frames (similarly for the longitudinal magnetization im- 
ages, which are not shown). This indicates that during 
this latency stage, the BEC remained in the unmagne- 
tized scalar phase (the \m z = 0) state). Any magnetiza- 



FIG. 2: In-situ images of ferromagnetic domains and domain 
walls, (a) The transverse magnetization density, measured 
for condensates at variable times Tn„jj (indicated on top), 
is shown with the magnetization orientation (<f> — arg(Fr)) 
indicated by hue, and magnitude (A = CI -Ft I) indicated by 
brightness. The maximum brightness, indicated by the color 
wheel at left, corresponds to full magnetization of the con- 
densate center. For the data at Ihold = 216 ms, the mag- 
netization density (b) and orientation (c) are shown sepa- 
rately, and regions with either a large (I) or several small 
(II) magnetized domains are indicated. The gray scale in (c) 
covers the range to 2n. Regions outside the condensate 
are indicated in black. Higher resolution picture available at 
http://physics.berkeley.edu/research/ultracold/pubs.html 

tion during this stage either was too low in magnitude 
or varied over too short a length scale to be discerned 
by our imaging. At later times, a non-zero transverse 
magnetization signal spontaneously developed, yielding 
a Larmor precession signal that varied both in ampli- 
tude and in phase across the condensate. This observa- 
tion indicates the spontaneous, spatially inhomogeneous 
breaking of 0(2) symmetry in the transverse plane in a 
direction given by the phase of Larmor precession. 

In Fig. 2, the derived transverse magnetization for 
samples at variable times Tjjo^ are presented, rendered 
as color images to portray both the magnetization orien- 
tation (as hue) and amplitude (as brightness). These 
spatial maps show ferromagnetic domains of variable 
size and orientation arising spontaneously following the 
quench. The landscape of domains includes small re- 
gions of homogeneous magnetization, with bands of un- 
magnetized gas, i.e. domain walls, separating regions of 
nearly opposite orientation. One also observes larger fer- 
romagnetic spin textures free of domain walls in which 
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the magnetization orientation varies smoothly over tens 
of microns. 

Ferromagnetic domains: The evolution toward 
ferromagnetism from the \m z = 0) state under the con- 
ditions q < \ c2\n is a form of dynamical instability in a 
spinor bec dEEHIil. This instabilit y 

is a conse- 
quence of coherent collisional spin mixing in which pairs 
of atoms in the \m z — 0) state collide inelastically to 
produce atoms in the \m z = ±1) states [!, 01 IS El ■ In- 
deed, in our work, measurements of populations in the 
different magnetic sublevels through Stern-Gerlach anal- 
ysis showed significant mixing of populations coincident 
with the onset of spontaneous Larmor precession. 

One may also regard this instability as the phase sepa- 
ration of a two-component condensate 0, ^g] j recalling 
that the \m z = 0) state represents an equal superposition 
of the 1 77i0 = ±1} eigenstates of any transverse spin oper- 
ator, Ftp — F x cos 4> + F y sin</>. Additionally, since C2 < 0, 
the ±1 eigenstates of any spin component are immis- 
cible [l^. Ferromagnetism thus arises by the spinodal 
decomposition of a binary gaseous mixture into neigh- 
boring regions of oppositely oriented transverse magne- 
tization. Applying results from the two-component case, 
this phase separation is dominated by an instability with 
a characteristic exponential timescale tfm = h/ *J2\c2\n 
with 77 being the total gas density. The wavevector of 
the dominant instability kpM = y/2m\c2\n/h, defines the 
typical size I = 7rkp M of single-component domains in 
the phase-separated fluid, and also the width b ~ k F \j of 
domain walls in which the two components still overlap 

MM- 

To relate these predictions to our data, we considered 
the density-weighted transverse magnetization correla- 
tion function 



(a) 



G T (Sp) = Re 



Sp (C(P)JHP))* (C(P + Sp)F T (p + 5p)) 



At zero range, Gt(0) measures the degree to which the 
condensate has evolved toward the ferromagnetic state. 
As shown in Fig. 3 (a), Gt(0) rises after the quench 
from a near-zero background level characteristic of the 
unmagnetized initial state of the condensate to a nearly 
constant value of Gt(0) — 0.5 during the ferromagnetic 
stage. This saturation value measures the area occu- 
pied by domain walls. Fitting Gt(0) at early times 
(^hold < 90 ms) to a rising exponential yields a time 
constant of r = 15(4) ms. This value is in agreement 
with the predicted tfm = 13.7(3) ms, calculated using 
the peak density averaged over the y direction to account 
for its reduced dimensionality. Over the same period of 
evolution, no significant longitudinal magnetization was 
observed, confirming the presence of purely transverse 
ferromagnetic domains. 

Spatial correlations in the transverse magnetization 
(Fig. 3(b)) are typified by a central region of positive 
correlations (near 5p — 0) and then several equally- 
sized regions of alternating negative and positive corre- 
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FIG. 3: Temporal and spatial evolution of ferromagnetism in 
a quenched spinor BEC. (a) The mean local squared trans- 
verse (Gt(0), squares) and longitudinal (Gl(0), circles) mag- 
netization, averaged over several experimental repetitions (er- 
ror bars indicate shot-to-shot rms fluctuations), is shown. 
Gt(0) rises exponentially with time constant r = 15(4) ms 
as determined by a fit to data for TvjqU < 90 ms (gray line), 
and saturates at Gt(0) ~ 0.5 at later times. No significant 
longitudinal magnetization is observed. The minimum val- 
ues for Gt(0) and Gt(0), evident in the data for earliest 
-^hold' renec t residual noise due to processing either 24 or 
2 frames, respectively, to obtain the magnetization density, 
(b) Maps of spatial correlations in the transverse magnetiza- 
tion (Gt(x, z)), at variable T no i(j show alternating regions of 
positive and negative correlations in the narrow x direction, 
indicating that phase separation occurs primarily along x via 
a few discrete dynamical instability modes. At each Tu Q yi, 
Gt(x, z) was determined for ten repetitions of the experiment 
and then averaged. 



lations displaced from one another in the narrow (x) di- 
mension of the condensate. We estimate a typical size 
for single-component domains as ~ 10 /im, twice the 
displacement at which the transverse spin-spin correla- 
tions changes sign, in good agreement with the predicted 
Ttkp M — 8.3(2) /jm. The presence of negative correla- 
tion regions supports the model of spin-conserving phase- 
separation discussed above. The preferential phase sep- 
aration in the narrow x direction rather than along z is 
due to the broader momentum distribution of the un- 
magnetized condensate in that direction due to the finite 
condensate size. Thus, upon quenching the system, a 
greater population of atoms is available to seed the faster- 
growing, shorter- wavelength instabilities with wavevector 
in the x direction. The presence of several alternations 
of positive and negative correlations further suggests that 
the phase separation occurs through a small number of 
discrete, unstable magnetization modes. 
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Topological defects: The spontaneous symmetry 
breaking observed in this work is one of many examples 
of symmetry breaking that occur in nature. Symmetry 
breaking is presumed to have occurred at thermal phase 
transitions in the early universe, giving rise to the spe- 
cific elementary particles and interactions observed in the 
present day. An important aspect of rapid spontaneous 
symmetry breaking, whether in the laboratory or of a cos- 
mological nature, is the creation of topological defects, a 
process described theoretically by Kibble and Zurek 
|23| . The types of topological defects that may be formed 
depend on the group structure of the ground-state man- 
ifold reached at the transition. 

Analogues of cosmological topological defect formation 
have been studied using liquid crystals [24| and superfluid 
helium [2^, |2(| |53 • In comparison with these previous 
experiments, the present investigation focuses on a sim- 
pler physical system, in which the time for thermal equili- 
bration and symmetry breaking is much longer than that 
needed to bring the system across the symmetry-breaking 
transition. Furthermore, the present work focuses on a 
quantum rather than a thermal phase transition. 

In our two-dimensional system, spin- vortices are topo- 
logical point defects about which the orientation of mag- 
netization has a 2im winding with n being a non-zero 
integer. Following a procedure for identifying such de- 
fects in our data (Supplementary Note 3), spin- vortex 
defects were observed with high confidence in about one 
third of all images containing significant ferromagnetism 
(^hold > 90 ms) , with some images containing up to four 
vortices. Data revealing one such spin-vortex are shown 
in Fig. 4. For this vortex, the central region of near-zero 
Larmor precession amplitude has a diameter of about 3 
/im, comparable to the spin healing length £ s = 2.4 /um 
defined earlier. All the vortices we observed were singly 
quantized, with no apparent preferred direction of circu- 
lation. 

Many types of vortices that can occur in a gas with 
a multi-component order parameter can be distinguished 
by the composition of their cores 28] . Based on our mea- 
surements of the transverse and longitudinal magnetiza- 
tion at the vortex core, the spin-vortices seen in our ex- 
periment appear to have unmagnetized filled cores. This 
observation supports their characterization as "polar- 
core" spin- vortices (denoted as (1,0,-1) vortices in Ref. 
|28|h for which the superfluid order parameter is a super- 
position of atoms in the \m z = 1} state rotating with one 
quantum of circulation, atoms in the \m z = —1) state 
rotating with one quantum of circulation in the oppo- 
site sense, and non-rotating atoms in the unmagnetized 
\m z — 0} state, which also fill the vortex core. Such a vor- 
tex is thus characterized by zero net mass circulation and 
a spin current with one quantum of circulation. The ori- 
gin for such spin currents is presumably the spinodal de- 
composition by which ferromagnetism emerges from the 
unmagnetized cloud. The alternate identification of these 
vortices as merons [2^, [3(3 is ruled out by the absence of 
a longitudinal magnetization signal at the vortex core. 
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FIG. 4: In-situ detection of a polar core spin-vortex. Spa- 
tial maps of the transverse magnetization (a) magnitude 
(A = CI^t|) and (b) orientation (<f> — arg(.FT)) are shown 
for a sample imaged at Tjj jj = 150 ms. Data from a portion 
of the condensate, indicated by boxes, are magnified, showing 
the transverse magnetization (c) magnitude and (d) orienta- 
tion and also (e) the magnitude of the longitudinal magnetiza- 
tion (Al = £\F Z \). The phase along a closed path (indicated 
in red), surrounding a region of near-zero transverse magne- 
tization and satisfying the criteria stated in the text, shows 
a net winding of 2-k, revealing the presence of a spin- vortex 
defect with its core within the closed path. The core shows 
no significant longitudinal magnetization, allowing the iden- 
tification of the defect as a polar core spin- vortex. The gray 
scale for images (a,c,e) is the same. The field of view is 31 x 
320 fim for images (a,b), and 26 x 26 /im for images (c-e). 

In contrast with their topological nature in some other 
magnetic systems, domain walls in a F = 1 ferromag- 
net are not topologically stable; rather, they may decay 
by the formation of spin vortex- antivortex pairs. As dis- 
cussed above, we observe domain walls that form at the 
onset of visible ferromagnetism in the gas, permeate the 
inhomogeneous clouds, and persist for all times there- 
after. The persistence of domain walls in our system is a 
subject for future investigation. 

Of keen interest for future studies is the distinction 
between thermal and quantum fluctuations as seeds for 
the unstable symmetry-breaking process. In this work, 
the forced depletion of atoms not in the \m z = 0) state 
suggests that ferromagnetism formed purely by the am- 
plification of quantum fluctuations, i.e. shot noise, a sug- 
gestion that warrants experimental justification. Time- 
resolved experiments in which one varies the rate at 
which the system is swept into the ferromagnetic state 
may also uncover universal temporal dynamics that typ- 
ify this and other quantum phase transitions. Our 
method of nondestructive magnetization imaging also al- 
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lows the behavior of a single quenched condensate to be 
studied in detail, permitting an examination of the dy- 
namical evolution of domain walls and spin-vortices. 

We thank E. Mueller, J. Moore, and A. Vishwanath 
for helpful comments, J. Guzman for experimental as- 
sistance, and the NSF and David and Lucile Packard 
Foundation for financial support. S.R.L. acknowledges 
support from the NSERC. 

Supplementary Note 1: Experimental sequence 
Optically-trapped BECs used for this work were obtained 
by loading around 10 8 atoms in the \F = l,m z = — 1} 
state with temperature of 2.5 [iK into an optical dipolc 
trap (ODT). The ODT was formed by a single focus laser 
beam with a wavelength of 825 nm that is linearly polar- 
ized to ensure that all components experience the same 
trap potential. The weak quadratic dependance of the 
trap depth on nip is negligible. After loading the atoms, 
all were placed in the \m z = 0) magnetic sublevel using 
rf rapid adiabatic passage followed by application of a 
transient magnetic field gradient of 4 G/cm. Holding the 
magnetic field at a magnitude of 2 G, the optical trapping 
beam was decreased in power over 400 ms, lowering the 
trap depth to fc^x 350 nK and the temperature of the 
gas to ~ 40 nK, well below the Bose-Einstein condensa- 
tion temperature. This yielded nearly-pure condensates 
of 2.1(1) x 10 6 atoms. 

After forming the condensate, the magnetic field was 
ramped within 5 ms to a value of 50 mG oriented in the 
z direction. The field was held at this setting for a vari- 
able time before the imaging sequence. Field gradients 
along the x and z directions were nulled to less than 0.2 
mG/cm. A Stern-Gerlach analysis of populations in each 
of the magnetic sublevels was applied to establish that 



the field ramp was sufficiently slow as not to alter the 
spin state of the condensate. Such an analysis bounded 
the populations in each of the \m z = ±1) spin states to 
be less than 0.3% of the total population, both before 
and right after the field ramp. 

Supplementary Note 2: Magnetization imaging 

The probe light for the phase contrast images was de- 
rived from a diode laser detuned 5 = -200 MHz from the 
5S 1/2 (F = 1) 5P 1/2 (F' = 2) transition in 87 Rb. The 
imaging sequence to determine the transverse magneti- 
zation consisted of a sequence of 24 frames for which the 
probe was pulsed at a strobe frequency of around 10 kHz. 
The duration of each individual pulse was around 1 /j,s. 

Supplementary Note 3: Image Analysis The de- 
termination of spin vortices in the magnetization images 
relied on the identification of a core of "dark" pixels, 
consistent with zero Larmor precession amplitude, which 
were surrounded by "bright" pixels with a finite ampli- 
tude and well-defined phase of Larmor precession. We 
identify vortices based on two criteria: (1) that there is 
an island of dark pixels, a candidate for the vortex core, 
which is surrounded entirely by bright pixels and that is 
at least two bright pixels away from nearby dark pixels, 
and (2) that the transverse magnetization traced along 
a closed loop through bright pixels surrounding the core 
have a non-zero net winding. The distinction between 
bright and dark pixels (at about one quarter of the max- 
imum Larmor precession amplitude) was chosen so as 
to eliminate false-positive vortex detections in simulated 
data taking into account the measured noise and resolu- 
tion of our imaging system, and using estimates for the 
width of contiguous domain walls. 
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